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ABSTRACT 

The Schweizer-Middleditch star, located behind the SN 1006 remnant and near its center in projec- 
tion, provides the opportunity to study cold, expanding ejecta within the SN 1006 shell through UV 
absorption. Especially notable is an extremely sharp red edge to the Si II 1260 A feature, which stems 
from the fastest moving ejecta on the far side of the SN 1006 shell — material that is just encountering 
the reverse shock. Comparing i/ST far- UV spectra obtained with COS in 2010 and with STIS in 1999, 
we have measured the change in this feature over the intervening 10.5-year baseline. We find that the 
sharp red edge of the Si II feature has shifted blueward by 0.19±0.05 A, which means that the material 
hitting the reverse shock in 2010 was moving slower by 44 ± 11 kms~^ than the material that was 
hitting it in 1999, a change corresponding to —4.2 ± 1.0 kms^^ yi'^^- This is the first observational 
confirmation of a long-predicted dynamic effect for a reverse shock: that the shock will work its way 
inward through expanding supernova ejecta and encounter ever slower material as it proceeds. We 
also find that the column density of shocked Si II (material that has passed through the reverse shock) 
has decreased by 7 ± 2% over the ten- year period. The decrease could indicate that in this direction 
the reverse shock has been ploughing through a dense clump of Si, leading to pressure and density 
transients. 

Subject headings: ISM: individual (SN 1006, SNR G327.6-f 14.6) — shock waves — supernovae: indi- 
vidual (SN 1006) — supernova remnants 



1. INTRODUCTION 

The idea of a reverse shock that can reheat expanding 
ejecta in a young supernova remnant (SNR) was i ntro- 
duced in a pair of coinpanio n papers by iGullI ()1973l ) and 
[Rosenberg & S cheued (Il973f). a i id wa s further articulated 
bv lMcKea (| 19 7"! and bv fGuH (|l975l ). It was rapidly in- 
corporated into the in terpretation of earl y X-ray data 
from young SNRs (e.g.. ICharles"eFaI1[l975h . and has be- 
come firmly entrenched into the canonical model of a 
young SNR. Briefly, the idea is that following a super- 
nova (SN) event, a blast wave moves outward through the 
surrounding interstellar or circumstellar material, which 
becomes swept up behind the expanding strong shock. 
Within this shell, a contact discontinuity forms between 
the expanding SN ejecta and the newly swept-up mate- 
rial, and an internal shock develops that works its way 
inward through the ejecta, reheating them and giving 
rise to strong X-ray and optical emission from material 
that is highly enriched in heavy elements. As the SNR 
evolves and the reverse shock works its way inward, it 
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should encounter ever-slower ejecta — a proposition that 
we test through the observations reported here. 

Despite this long pedigree, detailed observations of the 
development of the reverse shock over time in actual 
young SNRs have been sparse. The most detailed has 
been for Cas A: the Chandra images clearly showing the 
reverse shock and its encounter with ej ecta of different 
compositions (jHwang et al.l 120001 I2004D and the Hub- 
ble Space Telescope (HST) images showing the evolution 
of individ ual ejecta knots a s they encounte r the reverse 
shock (M orse et a l.l l2004: Fesen et al.l[20Tll ). In this pa- 
per, we report the results of an experiment designed to 
test a definite prediction of the reverse-shock model: that 
as the reverse shock moves inward through the freely- 
expanding ejecta, the velocity of material entering the 
reverse shock gradually decreases over time. The freely 
expanding ejecta are too cold to be readily detected in 
emission, but the velocity profile can be probed through 
absorption spectroscopy carried out in the UV, where 
ions commonly found in ejecta have strong resonance 
lines connecting to the ground state. One of very few 
objects where this technique has so far been employed is 
the remnant of the bright Type la supernova SN 1006 
(SNR G327.6+14.6), located at a distance of 2.2 kpc 
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and ^ 500 pc above the Galactic plane ([Winkler et alj 
I2003D . Located only 2.'8 arcmin south of the projected 
center of the the 15' radius SN 1006 shell, the UV-bright 
iSchweizer fc MiddleditchI (| 19801 . henceforth SM) star en- 
ables study of a "core sample" through the remnant Q 

The SM star is an OB subdwarf, probably not much 
farther away than SN 1006. It is relatively bright {V ^ 
16.74, B — V = —0.14 ), with low fore ground cxtinc- 
tion, E(B - V ) = 0.1 (|Schweizer fc Mid dlcditch 1980; 
iWu et al1 ll993^ . Its potential as a background "UV light- 
bulb" for probing the interior of SN 100 6 through absorp - 
tion spectroscopy was first exploited bv lWu et al.l ()1983f ) 
using the International Ultraviolet Explorer (lUE), with 
subsequen t stud i es ca .rried out by |F cscn ct al. (1988); 
IWu et all (|T99l IT997I ). and IBlair et al.i (11996). These 
spectra show strong Fe II lines at 2383 A and 2600 A that 
are broad and nearly symmetric, and in addition lines 
from Si II, Si III, and Si IV, all of which are also broad, 
but with p rofiles that are e ntirely at positive (redshifted) 
velocities (jWu et al.l II 9971 ). The most detailed UV spec- 
trum of the SM star was obtained from the Hubble Space 
Telescope [HST) in 1999 with the Space Telescope Imag- 
ing Sp ectrograph (STIS) and reported bv lHamilton et al] 
(|2007l hereafter HFB07), which also summarizes the pre- 
vious UV studies. 

The strongest feature in the 1150-1700 A band is broad 
absorption due to Si II 1260.4 A, all on the red side of 
the line, extending from a velocity of ~ -1 -2400 kms~^ to 
a sharp edge at -h7026 kms'^ (HFB07). iHamilton et~all 
()1997l |2007[) have modeled this profile as a combina- 
tion of both shocked and unshocked Si II ejecta within 
SN 1006: the highest redshifts result from ejecta that 
are still freely expanding on the far side of the shell, 
with the sharp edge due to material just before it en- 
ters the reverse shock. The reverse-shock model predicts 
that, over time, the shock encounters ever slower- moving 
ejecta, and hence the red edge should gradually move to 
shorter wavelength. We report here a second-epoch spec- 
trum of the SM star obtained with the Cosmic Origins 
Spectrograph (COS) on HSTiYidX, when compared with 
the 1999 STIS spectrum, demonstrates this phenomenon 
for the first time. 

2. OBSERVATIONS AND WAVELENGTH REGISTRATION 

For this study, we used as the first-epoch spectrum the 
far-ultraviolet echelle spectrum of the SM star obtained 
with STIS in July 1999 and described by HFB07. Briefly, 
the STIS observation used a 2" x 2" slit, the E140M grat- 
ing and the FUV-MAMA detector to cover wavelengths 
1150-^1700 A at a scale of 0.0122 A pixer\ The line 
spread function is non-Gaussian, with a narrow core and 
broad wings, and in the crucial region of the spectrum 
near 1300 A has a FWHM of 4.6 kms'^ 

For the second-epoch observation, we used the newly 
installed COS instrument, since it has a far higher 
throughput than STIS, though with spectral resolution 
that is almost four times lower. The observation was 
carried out on 2010 Jan 16 using the 2'.'5 diameter pri- 
mary science aperture and the G130M grating. A central 
wavelength setting of 1327 A and standard subsampling 

^ In addition to the SM star. lWinkler et al.l l l2005l ) have used two 
fainter background QSOs to probe additional Hnes of sight through 
SN 1006. 



(FP-POS = auto) were used to cover the range ~ 1170- 
1470 A, with a gap extending over 1317-1327 A, at a 
dispersion of 0.010 A pixeP^. The COS line spread 
function is also non-Gaussian, with FWHM 18 kms^^ 
near 1300 A fsee ISahnow eFalliMoi . for a summary of 
the COS on-orbit performance). Because of the higher 
throughput with COS, only a single orbit was required 
to obtain a spectrum with signal-to-noise similar to that 
obtained in eight orbits with STIS. Key aspects of the 
two observations are summarized in Table 1, and in Fig. 

1 we plot the spectra from both observations, binned to 
emphasize the broad features. The time separation be- 
tween the epochs is 10.48 yr. 

In order to measure the small shift in the wavelength 
of the absorption edge, it is crucial that the STIS and 
COS spectra be precisely registered in wavelength. In 
order to achieve more precise registration than can be 
done using the wavelength scale provided with pipeline- 
processed COS data, we can use the many narrow in- 
terstellar and stellar lines in the wavelength region near 
the edge. We have used a cross-correlation analysis (the 
IRAFH task fxcor), first to cross-correlate each of the 
four individual (subsampled) frames, made tiny wave- 
length shifts to best align these, and averaged them to 
optimize the total COS spectrum. We then used a similar 
cross-correlation analysis to determine the small offsets 
between the COS and STIS spectra on both the blue 
(1245 - 1264 A) and red (1292 - 1301.5 A) sides of the 
edge, and found that with the original wavelength scales, 
the narrow lines differ in wavelength by as much as 40 
mA, or 9.6 kms~^. The differences are systematic, so 
we were able to apply a simple linear transformation to 
the COS wavelength scale to achieve registration with 
the STIS one to within ±6 mA, or ±1.5 kms~^. In Fig. 

2 we plot the difference between measured wavelengths 
in the COS and STIS spectra for several narrow lines, 
both before and after the transforming the COS spec- 
trum. After registering the COS and STIS spectra, we 
find a clear shift in the red edge of the 1260 A absorption 
feature to shorter wavelength, which we have measured 
as described in the next section. 

3. MODEL FIT TO THE COS SPECTRUM 

As discussed bv IHamilton et al.l ()1997l I2007D . the un- 
usual shape of the broad, redshifted Si II 1260 A absorp- 
tion feature in SN 1006 can be understood as arising from 
a combination of unshocked and shocked Si II ejecta. The 
sharp red edge at ^ 7000 kms~^ marks the position of 
a reverse shock front, where unshocked Si II is "instan- 
taneously" decelerated and thermalized. The lifetime of 
the resulting shocked Si II against colhsional ionization 
is of the order of the age of the remnant, so it is not sur- 
prising that shocked Si II should survive to remain ob- 
servable. While separated three-dimensionally of course, 
both shocked and unshocked Si II lie along the line of 
sight. 

Fig. [3] shows the result of fitting the 1260 A feature 
observed with COS to a combination of unshocked and 
shocked Si II. The analysis is essentially the same as 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the AURA, Inc. under cooperative 
agreement with the National Science Foundation. 
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that of HFB07, to which the reader is referred for more 
detaiL Fig. [3] can be compared to Figure 4 of HFB07, 
which showed a similar fit to the same feature observed 
with STIS. 

In brief, each of the observed spectra were interpolated 
across narrow interstellar lines, and their continua fitted 
to a low-order polynomial. In the case of the STIS spec- 
tra of the SM star and the comparison star PG 0839-1-399 
analysed in 2007, the continuum was a quintic polyno- 
mial in log flux versus inverse wavelength fitted over un- 
absorbed regions from 1150 to 1700 A. The COS and 
STIS spectra are in generally excellent agreement (see 
f|5]), but there is a small but significant difference in the 
overall level of fiuxes. To compensate for this difference, 
we multiplied the COS spectrum by a factor linear in 
log flux versus inverse wavelength, varying from 1.03 at 
1190 A to 1.08 at 1290 A. 

Fig. [3] also shows the ratio of the interstellar-line- 
excised, continuum-corrected spectrum of the SM star 
observed with COS, to that of the comparison star PG 
0839-1-399 observed with STIS. The compari son star PG 
0839-f3 99 (iGreen. Schmidt, fc Liebertlll986D was identi- 
fled by iBlair et al.l (|1996l ) as having not only a similar 
temperature and gravity, but also similar photospheric 
abundances and low extinction. As remarked by HFB07, 
the spectrum of the comparison star is not as similar in 
detail to that of the SM star as one might have liked, but 
dividing the SM spectrum by the comparison one does 
help to clarify which features in the SM star spectrum 
are due to SN 1006. 

Three separate velocity parameters can be measured 
from the line profile shown in Fig. [3] the free-expansion 
velocity rs/t (shock radius Vs divided by age t) of un- 
shocked Si II at the reverse shock front, and the mean 
velocity v and one-dimensional velocity dispersion a of 
shocked Si II. Not surprisingly, in view of the good agree- 
ment between the STIS and COS spectra, and the similar 
analysis, the parameters measured from COS are quite 
similar to those measured from STIS, Table [2] 

The three velocities are related by shock jump condi- 
tions. As noted by IIFB07, the measured velocities point 
to essentially all of the shock energy going into heating 
ions, with little energy in heating electrons, consistent 
with other observational evidence (Laming et al.' '1996'; 
[Ghavamian et al. 2002; Vink et al. 2003; Vink 2005), 
and also with Uttle energy in accelerating particles (e.g. 
iKosenko et al.l[2011h . In the absence of electron heating 
or particle acceleration, the change Au in the velocity be- 
tween the unshocked and shocked Si II should equal the 
three-dimensional velocity dispersion v^cr of the shocked 
Si II, 

Aw= y-t; = V3cr . (1) 

If some of the shock energy went into electron heating or 
particle acceleration, then the velocity dispersion would 
be correspondingly smaller. If all three velocities are 
treated as free parameters, the observations mildly pre- 
fer the opposite, an ion velocity dispersion slightly (Icr) 
larger than predicted by equation In the present 

paper we choose to constrain the fit. Fig. |31 such that 
equation (jlj is satisfied, whereas HFB07 chose to mea- 
sure all three velocities separately from the STIS spec- 
trum. IIFB07 estimated a 3(t upper limit of 0.26 of the 



shock energy going into forms other than ion heating. 

Fits constrained to satisfy equation ^ prefer slightly 
less shocked Si II absorption at its peak near 
5,000 km s~^ . A simple way to accomodate the preference 
is to fit the blue edge of the profile up to 5,100 kms~^ in- 
stead of the 5,600 kms~^ adopted by IIFB07. The values 
quoted in Table [5] are best fits over the reduced range. 
The quoted uncertainties are la statistical uncertainties 
about the best fit, subject to variations not constrained 
by equation ([T]). 

The observed column density of shocked Si II mea- 
sured from the line p rofile in Fig. [3] is (compare eq. 8 of 
IHamilton etaIlll997D 

^Siii = "l^f \/2ti utq — 7.4 ± 0.3 X 10^^ cm~^ (measured), 
7re^ / A 

(2) 

wher e A = 1260 A is the wavelength of the line, / = 
1.007 fMo rtonlll99lD is its oscillator strength, tq is the 
optical depth at line center of the fitted Gaussian profile 
of shocked Si II, and the factor -s/^tt comes from integrat- 
ing over the profile. The uncertainty in the estimate ([5]) is 
the statistical uncertainty arising from photon shot noise, 
and does not include "systematic" uncertainty arising 
from modeling choices, which include placement of the 
continuum, and the range of velocities over which the 
fitting is done. The systematic uncertainty is somewhat 
larger. Notably, if we extend the range over which we fit 
the blue edge of the shocked Si II profile to 5,600 kms~^ 
(as in IIFB07) instead of 5,100 kms^^, then the mea- 
sured column density increases to 8.0 x lO^^cm"^. The 
measured value ^ is listed in Table [2l along with the 
value that HFB07 would have measured from STIS but 
did not report. 

The measured column density ^ of shocked Si II is 
consistent with that predicted by steady ionization of 
Si II (compare eq. 13 of .Hamilton et al.i,1997, ). 

,rshk nsilIVs 
^^SiH — -j — -, — r 

= 6±2x lO^^cm'^ ( "iT" ) '(P^s<iicted), (3) 

where (crw)siii is the collisional ionization rate (see be- 
low), and the quoted uncertainty is a nominal estimate of 
the uncertainty in that rate. The electron to Si II ratio of 
rie/nsiii = 1.5/0.6 — 2.5 in equation (jSj is a fiducial value 
that assumes that each Si ion contributes on average 1.5 
electrons during its steady state ionization, and that the 
Si is in a silicon-rich region where it comprises 60% of 
the ion abundance by number, the remaining 40% being 
mostly S, consisten t with SNIa nucleosynthesis models 

(Hofii ch et al1[l998hFI 

As discussed bv IHamilton et ahl (|1997| ) , it is easy to re- 
duce the theoretical value ^ of shocked Si II, by diluting 
the Si with other elements (i.e. by increasing ng/nsiii), 
or by truncating the Si downstream, but it is harder to 
increase the theoretical value, to approach the slightly 
greater observed value ([2])- Thus agreement between 
the measured and predicted column densities of shocked 
Si II appears to require that the absorbing shocked Si 

^ Many lines of evide nce suggest that SN 1006 was a Type la 
event Ce.g.. I Winkler et al. 2003 ). 
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arises from silicon-rich ejecta. This conclusion is consis- 
tent with the absence of Fe II absorption with the same 
profile as the shocked Si II. Although Fe II absorption is 
observed, its profile is dramatically different from that of 
the shocked Si II, indicating th at little iron is mixed in 
with the shocked sili c on (Hamilton et al.l 119971 ) . 

In iHamilton et al.l (|1997t ) there seemed to be a pos- 
sible discrepancy between the observed and predicted 
values of the column density of shocked Si II. That 
discrepancy has now disappeared, in part thanks to 
a reduction in the collisional ionization cross-section. 
The collisional ionization rate used in equation ([3|) is 
(cf )siii = 4.4 X 10~^ cm^s~^, obtained by in tegrat- 
ing cross-sections from iClark fc Abdallalil (|2003[ ) over a 
Maxwellian distribution of colliding electrons at temper- 
ature 95 eV, the temperature reached by electrons as a 
result of Coulomb collisions with Si II ions over the colli- 
sional ionizati on timescale of Si II . The c ollisional ioniza- 
tion rate from IClark fc Abdallahl (j2003l ) is smaller than 
the rat e 6.1 x lO'^cm ' ^s"^ from iLennon et all (| 19881 ) 
used bv lHamilton et al] ()1997D . The newer rate, based on 
theoretical computations, should be more reliable than 
the older rate, based on isoelectronic scaling. 

A concern is that the collisional ionization time ^siii of 
Si II is comparable to the age t of the remnant (specif- 
ically, tsiii/t ~ 0.8, given the pre-shock Si II density 
from the observed optical depth of freely expanding un- 
shocked Si II at the shock front, an estimated electron- 
to-Si II density rie/nsm ~ 2.5, and the theoreti cal col- 
lisional ionization rate (cf. IHamilton et allllQQTt )!. The 
collisional ionization time is long enough that the shocked 
Si II absorption should be a superposition of components 
shocked at various times with various velocities. It is 
somewhat surprising that the profile should fit well to a 
single component satisfying the shock jump conditions. 
A possible explanation of the coincidence is that one- 
dimensional hydrodynamic models predict rather flat ve- 
locity and temperature profiles of shocked ejecta (e.g. 
iDwarkadas fc Chevalieilll998 ' Koscn ko et al.l l201lV An- 
other possible explanation, suggested by the faster than 
expected change in the optical depth of shocked Si II, 
fj5l is that the absorption is being produced by a lump of 
silicon-rich ejecta whose average density was higher than 
that observed today, and which has been shocked over a 
timescale shorter than the age of the remnant. 

4. EVOLUTION OF THE SHOCK FRONT FROM STIS TO 

COS 

As the reverse shock propagates into the Si II ejecta, 
the free-expansion velocity rs/toi unshocked gas entering 
the shock at radius should decrease with time t as 

dir^^_v (4) 
at t 

where Vg is the velocity of the reverse shock relative to 
the ejecta. The shock velocity Vg is related to the change 
Av in the velocity between unshocked and shocked Si II 
by Vg = (4/3) Aw (regardless of electron heating). The 
deceleration Aw measured from the COS spectrum shown 
in Fig. [3]yields a shock velocity of = 2630±120kms^^, 
in good agreement with the STIS measurement of Vs — 
2680±120 kms-^ As described by HFB07, equation (g]) 
coupled with the measured shock velocity predicts that 
the free-expansion velocity at the shock should be de- 



creasing at the rate of — 2.7± 0.1 kms^"'^ yr"^. The pre- 
dicted change A(rs/t) in the free-expansion velocity of 
the shock over the 10.5 year interval between the STIS 
and COS observations is then 

A{rs/t) = -28 ± 1 kms"^ (predicted). (5) 

Fig. |4] compares the COS and STIS spectra in 
the vicinity of the putative reverse shock front near 
7000 kms^^. To permit direct comparison of the spec- 
tra, the COS spectrum has been convolved with the 
STIS line spread function (LSF), while the STIS spec- 
trum has been convolved with the COS LSF. The net 
resolution is essentially that of COS, since the COS res- 
olution (18 kms~^ FWHM) is over three times broader 
than STIS (4.6 kms^^ FWHM). The flux in the COS 
spectrum has been multiplied by a slowly-varying factor, 
about equal to 1.08 near 1290 A, to bring it to the same 
normalization as STIS, as described in ^ 

The spectra in Fig. |4] show that the shock front near 
7000 kms~^ has indeed moved during the interval be- 
tween the STIS and COS observations. One way to mea- 
sure the velocity change would be to shift the two spectra 
horizontally, but the plethora of stellar lines precludes a 
precise measurement by this naive approach. 

We follow instead a more rigorous approach that mod- 
els the physics of the situation and that is unbiased by 
the presence of stellar lines. Fig. [3] shows that, at the 
shock front, unshocked gas is attenuating the stellar flux 
by a factor of 0.4. If this is correct, and the shock front 
has shifted by A(r,/t) from STIS to COS, then the STIS 
spectrum should equal the COS spectrum multiplied by 
an inverted top hat 1 — 0.4 — 0.6 deep and A{rs/t) wide. 
We can think of the top hat as representing material that 
was absorbing in 1999 that has been "eaten away" by the 
reverse shock since. Observationally, this absorbing top 
hat will be convolved with the overall LSF of the spec- 
trum. 

Fig. |4] shows the best-fit COS spectrum adjusted by 
this procedure, with the convolved, absorbing top hat 
shown as a dotted line. The fit was obtained by minimiz- 
ing the between the STIS and adjusted COS spectra, 
the noise being taken to be photon shot noise. The cor- 
relations between spectral bins introduced by convolving 
the spectra were taken into account. 

The best fitting velocity shift between STIS (1999) and 
COS (2010) spectra is 

A(r,/t) = -44± 11 kms-\ (6) 
or — 4.2 ± 1.0 kms~^ yr~^ (measured). 

The measured 1 cr uncertainty of 11 kms^^ is about 60% 
of a COS resolution element. There can be little doubt 
that the velocity shift is real, as the measured value dif- 
fers from zero by 4 a. The measured velocity shift ([6]) is 
reasonably consistent with the predicted shift ([5]) , differ- 
ing by 16± 11 kms^^, or 1.5 cr. The value of 6982 kms~^ 
for rs/t given in Table [2] is just the STIS value (which is 
more precisely measured) reduced by 44 kms~^. 

It should be remar ked that we have a ssumed in th is pa- 
per, as previously in IHamilton et al.] (|1997l I2007D . that 
the shock front is perpendicular to the line of sight, a 
plausible assumption given that the SM star is close to 
the projected center of the SN 1006 remnant. If the shock 



Reverse Shock Evolution in SN 10061006 



5 



front normal were inclined by an angle 9 to the line of 
sight, but the ejecta are moving directly away from us, 
then the main effect would be to increase the observed 
rate of change of the free-expansion velocity by a factor 
1/ cosO. If a future observation should measure a change 
faster than that predicted by equation ([5]) by a statis- 
tically significant margin, then one natural explanation 
would be that the shock front is not perpendicular to the 
line of sight. 

5. APPARENT DECREASE IN SHOCKED Si II 
ABSORPTION 

An unanticipated change between the two spectra is 
that the optical depth in shocked Si II appears to have 
decreased over the 10.5 year interval between STIS and 
COS observations. The change is illustrated in Fig. [5l 
which compares the COS (2010) and STIS (1999) spectra 
over the full range of the broad, redshifted Si II absorp- 
tion feature. The ratio of COS to STIS fiuxes indicates 
a small but significant flux excess, with a profile similar 
to that of the shocked Si II illustrated in Fig. [S] As dis- 
cussed in ^iSl the COS flux has first been multiplied by 
a linear function varying from 1.03 at 1190 A to 1.08 at 
1290 A in order to bring its overall flux to the same level 
as STIS. 

To measure the apparent change in absorption, wc first 
adjusted the COS spectrum for the 44 kms^^ shift in the 
shock front as described in 21 so that the two should be 
comparable. We then assumed that the COS:STIS ratio 
has the same profile as shocked Si II but with unknown 
amplitude, and did a one-parameter least-squares fit to 
determine that amplitude. This procedure is more pre- 
cise than simply taking the ratio of shocked column den- 
sities in Table [21 whose values are model-dependent. We 
find that the column density of shocked Si II measured 
from the COS spectrum is 0.93 ± 0.02 times that in the 
STIS spectrum from ten years earlier. The statistically 
significant 7 ±2% change is larger than the ~ 1% change 
one might simplistically have expected over 10 years of 
observation of a 1000 year old supernova remnant. The 
2% uncertainty includes only statistical uncertainty, and 
not any systematic uncertainty associated with our mod- 
eling of the continuum. The 7% change is comparable to 
the ~ 8% systematic adjustment we have made to the 
COS flux, and some systematic uncertainty attends our 
modeling of the relative sensitivities between COS and 
STIS. 

If the unexpectedly rapid decrease in optical depth 
of shocked Si II is real, then a possible explanation is 
that the shock has been ploughing through a clump of 
Si. Numerical simulations of shocks overrunning inter- 
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stellar clouds (jKlein et al.lll994[ ) indicate that clouds of 
sufficiently high density contrast can undergo significant 
pressure and density transients as the shock propagates 
through the cloud. The line profile of unshocked Si II 
illustrated in Figure |3| does indicate a density decreas- 
ing inward by a factor of 2 over an interval of about 
700 kms^^. This change is equivalent to about 2% in 
30 kms^^, which is several times smaller than required 
by the observed change in optical depth. The profile of 
unshocked Si II shown in Figure [3] does show possible 
smaller scale structure consistent with more rapid varia- 
tion, but the reality of such structure is confused by the 
many stellar lines. 

6. SUMMARY 

Using UV absorption-line spectra against the SM star 
taken with STIS in 1999 and with COS in 2010, we have 
probed changes in the absorbing supernova ejecta seen 
along a "core sample" through the SN 1006 remnant. 
At both epochs we see strong absorption due to Si II, 
all of it red-shifted and with a sharp red edge. Follow- 
ing |Hamilto^£tIal] (1997, 20Q7,), we have interpreted the 
profile of the strong 1260 A absorption feature as result- 
ing from both unshocked and shocked (by the reverse 
shock) Si II ejecta, with the sharp red edge to this fea- 
ture being due to freely expanding ejecta just before be- 
ing hit by the reverse shock. After carefully taking into 
account small systematic differences in both wavelength 
and fiux between measurements with the two instru- 
ments, we find that the red absorption edge has shifted 
by — 44 ± 11 kms"'^ over the 10.5 yr interval between 
observations, or — 4.2 ± 1.0 kms^^ yr^^- This measure- 
ment confirms the long-held expectation that in a young 
SNR, the reverse shock will encounter ever-slower ma- 
terial as it works its way inward through the supernova 
ejecta. We also find the unanticipated result that the 
optical depth of shocked Si II has decreased by 7 ± 2%, 
faster than might be expected to result from homologous 
expansion over 10 years in the SNR's 1000-yr history. 
Clumped Si ejecta provides a possible explanation both 
for the decrease in optical depth and for the absence of 
corresponding blue-shifted Si II absorption. 



We appreciate the careful reading and thoughtful sug- 
gestions by the (anonymous) referee, which have helped 
clarify this paper. Support for this research has been 
provided by NASA Grant HST-GO-11659 from the Space 
Telescope Science Institute, which is operated by AURA, 
Inc. under NASA contract 5-26555. P.F.W. acknowl- 
edges additional support from NSF grant AST-0908566. 



Green, R. F., Schmidt, M., & Liebert, J. 1986, ApJS, 61, 305 
Gull, S. F. 1973, MNRAS, 161, 47 
— . 1975, MNRAS, 171, 263 

Hamilton, A. J. S., Fesen, R. A., & Blair, W. P. 2007, MNRAS, 

381, 771 (HFB07) 
Hamilton, A. J. S., Fesen, R. A., Wu, C.-C, Crenshaw, D. M., & 

Sarazin, C. L. 1997, ApJ, 481, 838 
Hoflich, P., Wheeler, J. C, & Thielemann, F. K. 1998, ApJ, 495, 

617 

Hwang, U., Holt, S. S., & Petre, R. 2000, ApJ, 537, L119 
Hwang, U., et al. 2004, ApJ, 615, L117 

Klein, R. I., McKee, C. F., & Colella, P. 1994, ApJ, 420, 213 



6 



Winkler, et al. 



TABLE 1 

Observations of the Schweizer-Middleditch star from HST 



Instrument 


Obs. Date 


Grating 


FWHM 
( kms-i) 


Exposure 
(s) 


HST Program 
ID 


STIS 
COS 


1999 Jul 24-25 
2010 Jan 16 


E140M 
G130M 


4.6 
18. 


23,186 
1,756 


7349 
11659 





TABLE 2 

Measured parameters 



Parameter 


COS (this paper) 


STIS (HFB07) 


rs/t 


6982 ± 11 kms-i 


7026±3(rel.)±10(abs.) kms"! 


V 


5030 ± 70 kms-i 


5160 ± 70 kms-i 


a 


1140 ±50 kms-i 


1160 ± 50 kms-i 


Vs 


2630 ± 120 kms-l 


2680 ± 120 kms-i 




7.4 ±0.3 X 10"cm"2 


8.6 ±0.3 X 10i*cm-2 



8x10" 








1150 1200 1250 1300 1350 

Wavelength (A) 



1400 



1450 



1500 



Fig. 1. — The complete COS and STIS FUV spectra of the SM star. The absorption feature at 1217 A has contributions from geocoronal 
Lya emission, stellar and interstellar Lya absorption, redshifted Si II 1190.4 A and 1193.3 A absorption, and redshifted Si III 1206.5 A 
absorption. The most prominent "clean" feature is the broad, redshifted Si II 1260.4 A line, due to a combination of shocked and unshocked 
Si II ejecta. The COS spectrum has been displaced vertically for clarity. 
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Fig. 2. — The difference between measured wavelengths of several narrow lines, on both blue and red sides of the Si II absorption feature, 
is plotted as a function of wavelength. A simple linear transformation in the COS wavelength scale brings the two spectra into registration, 
with an RMS uncertainty of ±1.5 kms"-"^. 
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Fig. 3. — COS absorption spectrum of SN 1006 around the redshifted Si II 1260 A feature, showing the best-fit Gaussian profile of shocked 
Si II, and the residual unshocked Si II. The spectrum is the ratio of the interstellar-line-excised, continuum-corrected COS spectrum of the 
SM star to the STIS spectrum of the comparision star PG 0839-1-399, with each smoothed to a resolution of 80 kms~^ FWHM before their 
ratio was taken. The upper spectrum shows the interstellar-line-excised, continuum-corrected comparison stellar spectrum of PG 0839-1-399 
at the same resolution, shifted by -|-32 kms^^ to mesh its stellar lines with those of the SM star, and offset vertically to separate it from 
the SN 1006 spectrum. 
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Fig. 4. — Comparison of COS (blue) and STIS (red) spectra of the SM star behind SN 1006. To permit direct comparison, each spectrum 
has been convolved with the line spread functon of the other. Also shown is the best-fit adjusted COS spectrum (black), in which the COS 
flux has been attenuated by a factor of 0.4 over a slice of velocity 44 kms"'^ wide starting at the 7026 kms"'^ velocity of the shock front 
observed with STIS. The dotted line shows the factor by which the COS spectrum has been multiplied, an inverted top hat convolved with 
the overall LSF. The residual difference between the COS and STIS spectra is shown. 
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Fig. 5. — Comparison of COS (blue) and STIS (red) spectra of the SM star behind SN1006. Each spectrum has been convolved with 
the line spread function of the other. The lower panel shows the ratio of the COS to STIS fluxes. The smooth (cyan) line is the ratio 
smoothed with a gaussian to 320 kms~^ FWHM. The ratio shows a small but significant excess with a profile similar to that of the 
shocked Si II illustrated in Figurc|3] indicating a decrease in the optical depth of shocked Si II in COS compared to STIS. 



